Quantum optics of a Bose-Einstein condensate coupled to a quantized light field 
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We consider the interaction between a Bose-Einstein condensate and a single-mode quantized light 
field in the presence of a strong far off-resonant pump laser. The dynamics is characterized by an 
exponential instability, hence the system acts as an atom-photon parametric amplifier. Triggered by a 
small injected probe field, or simply by quantum noise, entangled atom-photon pairs are created which 
exhibit non-classical correlations similar to those seen between photons in the optical parametric 
amplifier. In addition, the quantum statistics of the matter and light fields depend strongly on the 
initial state which triggers the amplifier. Thus by preparing different initial states of the light field, 
one can generate matter waves in a variety of quantum states, demonstrating optical control over the 
quantum statistics of matter waves. 

PACS numbers: 03.75.-b,42.50.-p,42.50.Ct,42.50.Dv 



I. INTRODUCTION 



In many ways, recently developed Bose Einstein con- 
densates (BEC) of trapped alkali atomic vapors |l|,|[ are 
the atomic analog of the optical laser. In fact, with the 
addition of an output coupler, they are frequently re- 
ferred to as 'atom lasers' [|1. Despite many interesting 
and important differences, the chief similarity behind the 
analogy is that both optical lasers and atomic EEC's in- 
volve large numbers of identical bosons occupying a sin- 
gle quantum state. As a result, the physics of lasers and 
BEC involves stimulated processes, which due to Bose 
enhancement often completely dominate the spontaneous 
processes which play central roles in the non-degenerate 
regime. 

Just as the discovery of the laser led to the develop- 
ment of nonlinear optics, so too has the advent of BEC 
led to remarkable experimental successes in the once the- 
oretical field of nonlinear atom optics |Q^. Nonlinear 
optics typically involves the study of multi-wave mixing, 
epitomized by phenomena such as parametric down con- 
version and phase conjugation. Due to the presence of 
collisions, the evolution of the atomic field is also nonlin- 
ear, and multi-wave mixing has been predicted ||9Hl4| and 
observed in multi-component condensates Pq] , as well as 
in scalar condensates [|6|. 

At the root of most optical phenomena is the dynam- 
ical interaction between optical and atomic fields. Un- 
der certain circumstances, one can formally eliminate the 
dynamics of the atomic field, resulting in effective inter- 
actions between light waves. Under a different set of 
conditions, one can eliminate the electromagnetic field 
dynamics, resulting in effective atom-atom interactions. 
These are the regimes of nonlinear optics and nonlinear 
atom optics, respectively. These regimes, therefore, rep- 
resent limiting cases, where either the atomic or optical 
field is not dynamically independent, and instead follows 
the other field in some adiabatic manner which allows for 
its effective elimination. 



Outside of these two regimes the atomic and optical 
fields are dynamically independent, and neither field is 
readily eliminated. In this paper we investigate the dy- 
namics of coupled quantum degenerate atomic and op- 
tical fields in this intermediate regime. In particular, 
we investigate a system which is analogous to the non- 
degenerate optical parametric amplifier (OPA) |17| , p^ . 
However, whereas the OPA involves the creation of cor- 
related photon pairs, this system involves the generation 
of correlated atom-photon pairs. The purpose of this pa- 
per is to develop a detailed theory for the interaction of 
quantized atomic and optical fields, with emphasis on the 
manipulation and control of their quantum statistics and 
the generation of quantum correlations and entanglement 
between matter and light waves. 

The specific system we consider consists of a Bose- 
Einstein condensate driven by a strong far off-resonant 
pump laser which interacts with a single mode of an op- 
tical ring cavity counterpropagating with respect to the 
pump. The strong pump laser is treated in the usual 
manner as a classical, undepleted light field and further- 
more it is assumed to be detuned far enough away from 
resonance that spontaneous emission may be safely ne- 
glected. The cavity field, henceforth referred to as the 
'probe', is assumed to be weak relative to the pump, 
and is treated fully quantum mechanically as a dynamical 
variable. It is the dynamical interplay between this probe 
field and the atomic field which is the subject of interest. 
The pump serves as a sort of catalyst, inducing a strong 
atomic dipole moment, thus significantly enhancing the 
atom-probe interaction. 

Assuming that the probe field begins in or near the 
vacuum state, and the atomic field consists initially of 
a trapped BEC, the initial dynamics is dominated by 
a single process: the absorption of a pump photon by a 
condensate atom followed by the emission of a probe pho- 
ton. We remark that in the far off-resonant configuration, 
this is a two-photon virtual transition in which the ex- 
cited atomic state population remains negligible. Due to 



atomic recoil, the absorption/emission process transfers 
the atom from the condensate ground state to a new state 
that is shifted in momentum space by the two-photon 
recoil. This new state constitutes a second condensate 
component, which can be considered as a momentum side 
mode to the original condensate. [| As this side mode is 
populated, it begins to interfere with the original conden- 
sate, resulting in fringes [ pi| . These fringes are seen by 
the pump and probe fields as a spatial density grating, 
which then enhances the photon scattering process. 

This interplay between interference fringes and scatter- 
ing can act as a positive feedback mechanism, in which 
case the system is unstable, and is characterized by ex- 
ponential growth. Any small signal, including quantum 
noise, will be sufhcient to trigger the instability, result- 
ing in the generation of exponentially growing side mode 
and probe fields. Of course this exponential growth is 
eventually reversed by high intensity effects, so that the 
long-time dynamics is characterized by large-amplitude 
nonlinear oscillations. 

At the present time, we focus on the small-signal 
regime, characterized by exponentially growing fields. In 
this regime, we demonstrate that the quantum state of 
the probe and side mode fields depends strongly on the 
initial conditions, so that, e.g., by injecting a small co- 
herent light field into the probe, one can create an en- 
tirely different quantum state than that generated from 
the amplification of quantum vacuum fluctuations. The 
differences are manifested in both the quantum statistics 
of the individual field modes, as well as in non-classical 
correlations and entanglement between them. 

This rest of this paper is organized as follows. Section 
II gives the background and relates the current theory 
to previous works in a variety of fields. Section III out- 
lines the basic model for a quantized many-body atomic 
field interacting with a strong classical pump laser and a 
quantized optical cavity mode. In Sec. IV, coupled-mode 
equations are developed for the condensate and its mo- 
mentum side modes. These equations are then linearized 
in Sec. V, resulting in a three-mode model which is ex- 
actly solvable. Section VI then discusses the exponential 
instability, with emphasis on the effects of collisions. In 
Sec. VII the quantum statistics of the atomic and electric 



^This new state may or may not be in the same internal 
ground state as the original condensate, depending on the 
polarizations of the pump and probe photons. If the magnetic 
sublevels are different, then the transition would be termed a 
Raman transition, if the sublevels are the same, than it may 
be thought of as Rayleigh scattering or two-photon Bragg 
scattering M.pof. As the states are already distinguished by 
their center-of-mass momentum states, to further distinguish 
them by an additional quantum number would add nothing. 
Our model deals specifically with the Bragg scheme, however, 
with only minimal modifications it could be applied to the 
Raman scheme as well. 



fields are investigated, and the extent to which they can 
be manipulated is determined. In Sec. VIII atom-photon 
entanglement is discussed, including an examination of 
two-mode squeezing between atomic and optical fields 
Lastly, section IX is a discussion and conclusion, which 
includes estimates of the important physical parameters, 
as well as potential experimental obstacles. 



II. BACKGROUND 



The system we describe is in fact an extension into the 
ultracold regime of the theoretical work of Bonifacio and 
coworkers on the Collective Atomic Recoil Laser (CARL) 
1^,^. The original CARL theory treated the atomic 
center-of-mass motion classically, an approximation cer- 
tainly valid for hot atoms, but not sufficient to describe 
ultracold samples such as EEC's. Within this framework, 
the feedback mechanism which gives rise to the exponen- 
tial instability in the CARL was outlined using a slightly 
different, but complementary physical picture, where the 
classical atomic center-of-mass motion in the optical po- 
tential of the counterpropagating pump and probe fields 
is responsible for the grating formation. The theory was 
extended to the limit of zero temperature by assuming 
that all of the (classical) atoms begin from rest, leading 
to the discovery of the so-called 'CARL cubic equation', 
which gives the exponential growth rate of the instabil- 
ity in terms of the relevant system parameters. Out of 
a desire to better understand the quantum statistics of 
the probe field, an attempt at a quantum T = theory 
was made ||2J]. However, this attempt explicitly assumed 
that the wavefunctions of the individual atoms could be 
localized in both momentum and position space to an 
extent which violates the Heisenberg uncertainty princi- 
ple. Thus rather than being a true quantum theory, it 
still treated the atoms as following 'classical' trajectories, 
but now with small 'quantum' fluctuations included. 

Both the original classical CARL model, as well as the 
later 'quantum' model, fall within the ray-optics approx- 
imation for the atomic field. Clearly then, one would 
expect such models to break down as soon as the atomic 
dcBroglie wavelength becomes comparable to the period 
of the optical potential formed by the pump and probe 
fields. As the wavelength of the optical potential is twice 
the optical wavelength, this breakdown should occur near 
the atomic recoil temperature, which for typical alkali 
atoms is on the order of microKelvins. As subrecoil tem- 
perature atomic vapors are achieved routinely through 
a variety of cooling techniques, a theory which properly 
treats the quantum motion of the atoms is required if one 
desires to investigate the behavior of the CARL in this 
regime. 

With the ultimate intent of extending the CARL the- 
ory into the BEC regime, so that the unique coherence 
properties of condensates might be further understood 
and exploited by the interaction with dynamical light 



fields, a quantum model of the atomic motion was for- 
mulated |25|] , where it was confirmed that the ray-optics 
versions did indeed break down for temperatures of the 
order of the recoil temperature or below. In fact, at 
T = 0, a second threshold for the existence of the ex- 
ponential instability was discovered, occurring when the 
bunching process is overcome by matter-wave diffraction. 
For T > Tfi, however, it was shown that the previ- 
ous theories make indistinguishable predictions from the 
quantum theory. We remark that while in this work the 
atomic center of mass motion was treated quantum me- 
chanically, the light fields were still treated classically, 
hence predictions concerning the quantum statistics of 
either the atomic or optical fields could not be made. 

A full quantum model of both the atomic and optical 
fields was recently outlined in Eq], where the subjects 
of manipulating quantum statistics and atom-photon en- 
tanglement were first addressed. The present paper is a 
detailed elaboration and extension of that work, includ- 
ing significant new physics. For example, utilizing the 
familiar s-wave scattering approach of BEC theory, the 
effects of atom-atom collisions are incorporated into the 
CARL theory for the first time. Also, in an extension 
of the OPA analogy, the existence of two-mode squeez- 
ing is shown to occur between a condensate side mode 
and the probe optical field. The current approach also 
differs from earlier work in that the familiar spontaneous 
symmetry breaking technique is no longer applied to the 
condensate. Instead it is assumed that a condensate well 
below the critical temperature is better described by a 
number state than a coherent state, as recent work ap- 
pears to demonstrate |27|j2q| . 

The fully quantum model is similar in many ways to 
a system studied by Zeng and coworkers [p9| , in which 
the principle of manipulating the quantum statistics of 
a condensate by its interaction with a quantized light 
field was first proposed. This paper, however, does not 
recognize the existence of unstable (exponential) solu- 
tions nor the fact that the system can be triggered from 
quantum noise. We note that the unstable (exponen- 
tial) solutions, and the possibility to initiate them from 
quantum vacuum fluctuations are both crucial compo- 
nents of this present work. Lastly, we mention the con- 
nection to recent work on matter-wave amplification by 
Law and Bigelow [Q, which also explores the interac- 
tion between condensates and quantized light fields. In 
that work, however, the light field is assumed to be heav- 
ily damped, thus allowing for its dynamical elimination. 
As a result, only the properties of the atomic field are 
studied in detail. 

CARL theory, including the present version, is also 
closely related to the theory of Recoil Induced Reso- 
nances (RIR) pl[ , in which the effects of atomic recoil on 
the pump-probe spectroscopy of an atomic vapor is in- 
vestigated. This theory treats the atomic center-of-mass 
motion quantum mechanically. The probe field, however, 
is not typically treated as a dynamical variable. Hence, 
it does not include the effects of probe feedback, which 



are necessary for exponential behavior. A detailed com- 
parison of the RIR and CARL theories is given in [p2| . 



III. THE BASIC MODEL 

In this section we derive a fully quantized model of 
a gas of bosonic two-level atoms which interact with 
a strong, classical, undepleted pump laser and a weak, 
quantized optical ring cavity mode, both of which are as- 
sumed to be tuned far away from atomic resonances. As a 
result, single-photon transitions between atomic internal 
ground and excited states are highly non-resonant and 
the excited state population remains negligible. In this 
case, one can safely neglect the effects of spontaneous 
emission as well as the two-body dipole-dipole interac- 
tion. 

We must still, however, allow for two-photon virtual 
transitions in the which the atomic internal state remains 
unchanged, but due to recoil may result in a change in 
the atom's center-of-mass motion. For example, an atom 
which absorbs a pump photon and emits a probe photon 
experiences a recoil kick equal to the difference of the 
momenta of the two photons (which for nearly counter- 
propagating pump and probe beams is of the order of 
two optical momenta). These transitions, therefore, cou- 
ple different states of the atomic center-of-mass motion. 
Due to the quadratic dispersion relation of the atoms, 
these transitions will in general be non-resonant. For very 
cold atoms, the resultant detunings are typically on the 
order of the recoil frequency, i.e. much smaller than the 
natural linewidth of the atomic transition, ^a, whereas 
the one-photon transitions which we are neglecting have 
a detuning many orders of magnitude larger than 7^ . 

Our theory begins with the second-quantized Hamilto- 
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~r fiatom—pump i '^atom — atom: \^) 

where Hatom and Tiprobe give the free evolution of 
the atomic field and the probe mode respectively, 
Ti-atom^probe and Hatom-pump describc the dipole cou- 
pling between the atomic field and the probe mode and 
pump laser, respectively, and Hatom-atom contains the 
two-body s-wave scattering collisions between ground 
state atoms. 

The free atomic Hamiltonian is given by 



Ha 



= / d^r 



'^l^'^n-^^' + V9i^)]'^9i^) 



^lir){-—V' + nu;a + V,{r))Mr) 



, (2) 



where m is the atomic mass, to a is the atomic resonance 
frequency, !f'e(r) and ^g(r) are the atomic field opera- 
tors for excited and ground state atoms respectively, and 



Vg{r) and Ve(r) are their respective trap potentials. The 
atomic field operators obey the usual bosonic equal time 
commutation relations [!?j(r),!^|, (r')] = Sjji5'^{r — r'), 

and [<Fj{r),^j,{r')] = [ifj(r),#J,(r')] = 0, where j,/ = 

{e,g}- 

The free evolution of the probe mode is governed by 
the Hamiltonian 



n. 



probe 



= hckA^A, 



(3) 



where c is the speed of light, k is the magnitude of the 
probe wave number k, and A and A^ are the probe pho- 
ton annihilation and creation operators, satisfying the 
boson commutation relation [A,A^] = 1. The probe 
wavenumber k must satisfy the periodic boundary con- 
dition of the ring cavity, k — 2'Ki/ L, where the integer 
£ is the longitudinal mode index, and L is the length of 
the cavity. 

The atomic and probe fields interact in the dipole ap- 
proximation via the Hamiltonian 



T^a 



torn—probe 



ihgA / (frW^^{r)e'^"<l'g{r) 



H.c. 



(4) 



where g — d[ck/{2heoLS)Y^''^ is the atom-probe coupling 
constant. Here d is the magnitude of the atomic dipole 
moment, and S is the cross-sectional area of the probe 
mode in the vicinity of the atomic sample (where it is 
assumed to be approximately constant across the length 
of the atomic sample). 

In addition, the atoms are driven by a strong pump 
laser, which is treated classically and assumed to remain 
undepleted. The atom-pump interaction Hamiltonian is 
given in the dipole approximation by 



'i-atom—pump 



-.^e-'^"' f d^r<Pl{r)e'''"''<Pg{r) 



H.c, 



(5) 



where flo is the Rabi frequency of the pump laser, re- 
lated to the pump intensity /q by |rioP = 2d^Io/h eqc, 
LOo is the pump frequency, and fcg ~ wq/c is the pump 
wavenumber. The approximation indicates that we are 
neglecting the index of refraction inside the atomic gas, 
as we assume a very large detuning A = UQ — uja between 
the pump frequency and the atomic resonance frequency. 
Finally, the collision Hamiltonian is taken to be 



't-atom — atom — 



2Trh a 
m 



d=^r^t(j.)#t(r)^^(r)#,(r), (6) 



where a is the atomic s-wave scattering length. This cor- 
responds to the usual s-wave scattering approximation, 
and leads in the Hartree approximation to the standard 
Gross-Pitaevskii equation for the ground state wavefunc- 
tion (in the absence of the driving optical fields). 

We limit ourselves to the case where the pump laser 
is detuned far enough away from the atomic resonance 



that the excited state population remains negligible, a 
condition which requires that A ^ 7^. In this regime 
the atomic polarization adiabatically follows the ground 
state population, allowing the formal elimination of the 
excited state atomic field operator. We proceed by intro- 
ducing the operators ^^(r) = !^e(r)e''^"* and A' = Ae'-'^°\ 
which arc slowly varying relative to the optical driving 
frequency. The new excited state atomic field operator 
obeys then the Heisenberg equation of motion 



d ~ , , , 



zA#^(r) 



2 ^ 



f A\z.r 



^gir), (7) 



where we have dropped the kinetic energy and trap po- 
tential terms under the assumption that the lifetime of 
the excited atom, which is of the order 1/A, is so small 
that atomic center-of-mass motion may be safely ne- 
glected during this period. For the same reason, we are 
justified in neglecting collisions between excited atoms, 
or between excited and ground state atoms in the colli- 
sion Hamiltonian (ph. 

We now adiabatically solve for !f'g(r) by formally inte- 
grating Eq. (Q) under the assumption that !^g(r) varies 
on a time scale which is much longer than 1/A. This 
yields 
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1 

~ A 


^0 ikor 
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^g{r,0)e'^' 


+ ^e{ 


r,0)e^^*. 







(8) 



The third term on the r.h.s. of Eq. (pi) can be neglected 
for most considerations if we assume that there are no 
excited atoms at t = Q, so that this term acting on the 
initial state gives zero. The second term may also be ne- 
glected, as it is rapidly oscillating at frequency A, and 
thus its effect on the ground state field operator is neg- 
ligible when compared to that of the first term, which is 
non-rotating. ^ 

Dropping the unimportant terms, and then substitut- 
ing Eq. M) into the equation of motion for ^g(r), we 
arrive at the effective Heisenberg equation of motion for 
the ground state field operator 



dt 



^r) 



2m 

g|»ol 
2IAI 



Vg(r) Airha 



(ae'K-- -I- atg-.K.r^ 



^t(r)#g(r) 



We note that in much of the literature the second and third 
terms are simply ignored. We choose to keep them temporar- 
ily to demonstrate that the commutation relation for !f'e(r) 
is preserved (to order 1/A) by the procedure of adiabatic 
elimination. 



4A 



.9 

a' 



■^sW, 



(9) 



where K = k — ko is the recoil momentum kick the 
atom acquires from the two-photon transition, and we 
have introduced the new slowly-varying probe field oper- 
ator a = —i{flQA/\ilQ\\A\)A\ which stiU obeys the bo- 
son commutation relation [a,a^] — 1. Here, the second 
to last term is simply the optical potential formed from 
the counterpropagating pump and probe light fields, and 
the last term gives the spatially independent light shift 
potential, which can be thought of as cross-phase modu- 
lation between the atomic and optical fields. 

To complete our model, in addition to Eq. (O), we 
also require the equation of motion for the slowly varying 
probe field operator. By again substituting Eq. (||), we 
find that it obeys 



dt 



iS'a 



. g|»o| 
'2IAI 



j Sv^l{v)e-''^-^^g{v), 



(10) 



where 5' = lj{) — uj/is the detuning between the pump and 
probe fields. The probe frequency is given by w « ck, 
again assuming that the index of refraction inside the 
condensate is negligible. 



IV. COUPLED-MODE EQUATIONS 

We assume that the atomic field is initially in a Bose- 
Einstein condensate with mean number of condensed 
atoms N . Furthermore, we assume that N is very large, 
and that the condensate temperature is small compared 
to the critical temperature. These assumptions allow us 
to neglect the noncondensed fraction of the atomic field. 
Thus our model does not include any effects of conden- 
sate number fluctuations. 

We now introduce the atomic field operator which an- 
nihilates an atom in the condensate ground state 



CO = / d'v^l{v)^g{v), 



(11) 



where (/'o(r) = {Af^) satisfies the time-independent 
Gross-Pitaevskii equation 



_V2 - -i^ N\Mr)\^ + ^ (^o(r) = 0, 



2m 



% 



m 



(12) 



PL being the chemical potential. By differentiating Eq. 

(11) with respect to time, and inserting Eqs. (^ and 

(12) we find that the equation of motion for cq is 



1 ' f 



l^'' = -H^ + ^^^A«'"''° 
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''2IAI 



(fv~ 



^S(r)(iV|^o(r)p-'^j(r)^3(r))^g(r) 



(friflir) (ae'K-- + a) e-''^'') ^g{v). (13) 



From this equation we see that the effect of the optical 
fields is to couple the condensate mode to two side modes, 
whose wavefunctions are given by 



{r\ip±) = V3o(r)e 



ibiKr 



(14) 



In principle, the collision term in Eq. (n3) also cou- 
ples the condensate mode to various neighboring modes. 
However, to be consistent with the assumption of a pure 
condensate (T — 0) we assume that collisions alone do 
not populate any new atomic states. 

Defining the field operators for the (first-order) con- 
densate side modes as 



C±= / d'v{ip±\v)Wg{T), 



(15) 



allows us to reexpress the equation of motion for the con- 
densate mode field operator as 






d . . I ^J. 

"I 



m 



d'v^liv) N\^^{v)\'-m{v)^g{v) ^,(r 



2|A 



-{dc^ + a)c+), 



(16) 



where the operators Cj obey the bosonic commutation 
relations 



cj>c],] = {^jWj'); j,j' = {-,0,+}, 



(17) 



all other commutators being equal to zero. 

We note that the three states, \(po), and \(p±) are not 
mutually orthogonal, as their overlap integrals are given 

by 



(¥'Tl^±>-yrf'r|(^o(r)pe±'2^-'- 



(18) 



For most condensate sizes and trap configurations, how- 
ever, these integrals are many orders of magnitude 
smaller than unity. As a result, for 'typical' condensates, 
the orthogonality approximation 



iVjlV]') =%' 



(19) 



yields accurate results. The range of validity of this ap- 
proximation is discussed in Appendix A, where we briefly 
examine how the theory should be modified to properly 
take this non-orthogonality into account. In the follow- 
ing, however, we assume the validity of Eq. ([l^), so that 
the states \(po), and |(^±) can be considered as well de- 
fined and distinct modes of the atomic field. 

We now derive the Heisenberg equations for the mo- 
mentum side mode field operators, found by differentiat- 
ing Eq. ( p^ ) with respect to time and again inserting Eq. 
d), yielding 
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-I / d^r((j9_ 

m 



{^^Mv)(N\M^)\''-Wl{v)4',{v)Yg{v) 

(20) 
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2 A m 



where we have introduced four new field operators c±2 
and b±. 

The operators c±2, which have the definitions 



C±2 = / d?v{^±2\v)^g(v), 



(21) 



are the annihilation operators for the second-order side 
modes 



{^^±2) 



(/7o(r)e 



±i2KT 



(22) 



These modes will be optically coupled to third-order side 
modes, and so on so that a full theory of the nonlinear 
response of the system should include the entire mani- 
fold of side modes. In this paper, however, we focus on 
the linear regime, where only the first-order side modes 
contribute significantly. 

The operators h± have the definitions 



h±= / rf=^r(0±|r)#g(r), 



(23) 



where 



(r|0±) = (i^fc™,)-ie±*^"( ± zK • V(po(r)). (24) 

Here kc is the momentum width of the condensate state 
along K, and is roughly given by kc '^ 1/Wc, where Wc is 
the size of the condensate along K. The factor (Kkc)~^ 
is simply a normalization coefficient. 

To understand the physical meaning of the b± term in 
Eq. (20), consider what happens to a single atom after 



it is transferred into the state V'(r) = Vo(r) exp(— zK • r) 
at time t = 0. Under free evolution the wavepacket of 
the atom, which initially has the shape of the condensate 
ground state, will move with group velocity KK/m and 
spread at the velocity hkc/m. This evolution is described 
by the propagation equation 



V'(r,i) =exp[it(n/2TO)V^]i/'(r,0), 



(25) 



which for short enough times becomes 

{r\ip- 



V'(r, i) « ( 1 - it^) {r\ip-) +t-{K- V^o(r))e-^K- 



2m 



+ ^t^(vVo(r))e--•-. 



(26) 



The first term on the r.h.s. of Eq. ( pq ) gives a phase 
shift due to the kinetic energy of the atom, the second 
term contributes an infinitesimal translational shift, and 



the third term gives an infinitesimal amount of spread- 
ing. If we include the effects of the trap potential and 
collisions, this last term vanishes as all spreading effects 
are balanced by the trap potential for the ground state 
93o(r). From Eqs. (pq) and ( p4| ) we see that the state 
of the atom at time t can then be viewed as a coherent 
superposition of the state |(^_) and the state \<j>-). Thus 
the coupling to 6_ in Eq. (EQ) corresponds physically 
to translational nrotion of the side mode wavepacket at 
the recoil velocity Vr = TiK/m. Since the probability at 
time t that the atom is still in the ground state is simply 
the overlap between V(r,i) and V'('*, 0), it is clear that 
for times t <C Wc/vr this probability will be essentially 
unity, and the coupling to 6_ can be ignored. 



V. LINEARIZED THREE-MODE MODEL 

From Eq. (pol), we see that the first-order side modes 
are optically coupled to both the condensate mode and 
to second-order side modes. For times short enough that 
the condensate is not significantly depleted, the coupling 
back into the condensate is subject to Bose enhancement 
due to the presence of ~ iV identical bosons in this mode. 
The coupling to the second-order side mode, in contrast, 
is not enhanced. Hence for these time scales, the higher- 
order side modes are not expected to play a significant 
role. In addition, we consider only times t <^ Wc/vr, so 
that the translational coupling can be neglected. 

These arguments suggest developing an approach 
where the three atomic field operators cq, c_, and c+ play 
a predominant role. Therefore, we expand the atomic 
field operator as 

^gir) = (r|(^o)co + {r\^~)c- + {r\ip+)c+ + ^(r), (27) 

where the field operator 'ip{r) acts only on the orthogonal 
complement to the subspace spanned by the state vectors 
\(pa), |</5-), and \(p+)- As a result, ip{r) commutes with 
the creation operators for the three central modes. 

In the next step, we use Eq. ( |27| ) to expand the atomic 
polarization and collision terms in Eqs. (CT), (|l^), and 
(pO|), with the eventual goal of deriving a closed set of 
operator equations which fully describes the system dy- 
namics. At present, we are considering four dominant 
modes, the condensate and first-order side modes, as well 
as the optical probe mode. In the linear regime, however, 
we will see that the condensate mode can be dynamically 
eliminated, resulting in an effective three-mode model. 

In expanding the polarization and collision terms by 
means of Eq. (p7|), there are two principal considera- 
tions in determining which are the dominant terms. The 
first is Bose enhancement, which, in the regime of neg- 
ligible condensate depletion, strongly selects transitions 
involving the condensate mode. In order to estimate this 
effect, we assign a weight of ^/N for each occurrence of 
the operators cq and Cq in a given term. The second con- 
sideration is momentum conservation, which comes from 



the spatial integration in the polarization and collision 
terms. Integrals over slowly varying functions such as 
|(/3o(r)p, or |(po(r)|'* are 'momentum selected' and dom- 
inate over integrals of rapidly oscillating functions such 
as |iy9o(r)P exp(— iK • r). 

With this approach we find that the equation of motion 
for the probe field operator ( [lO| ) becomes 



dt 



iS'a 



. .9l»o| 
'2IAI 



cLco 



■%'^+ 



(28) 



Thus we see that the probe annihilation operator is cou- 
pled to the bilinear atomic field operators c_Co and CqC+. 
These operators correspond physically to interference 
fringes, i.e. a periodic modulation of the atomic density, 
which appear because the atoms are in a coherent super- 
position of the side mode and condensate states. Gain in 
the probe can thus be interpreted as Bragg scattering of 
the pump due to the presence of interference fringes. 

By inserting (p^) into Eq. (Hq) we further find that 
the equation of motion for cJcq is given to leading order 
in the collision and optical terms by 



dt 



CqCo 






2|A| 



Co 



CqC^ 



H.c. 



where 



Fo 



d'r\Mr)\' 



(29) 



(30) 



Similarly we find that the operators c_co and CqC+ obey 
the equations 



d 4 . _ . hK^ 
dt ^ 2m 



clcq + i- 



. AnhaFo 



, . g|»oPL 4^ 



and 



^4, 
Jt'^^' 



i 

2m 



4c+ 



. AnhaFo 






(cLco + cJc+)4co 

(31) 

(cLco + cJc+)cJco 

(32) 



We assume that all N atoms are initially in the con- 
densate mode, so that 



1 / A^ 
|^),.o = ^(4) |0), 



(33) 



cjco 



Nil + So), 



cLco 



iV(5_, 



and 



sts _ 



NS 



+ 1 



(34) 



(35) 



(36) 



where Sq, (5_ and (5_|- are therefore infinitesimal operators. 
In addition, we introduce a rescaled probe field operator 



a 

N' 



(37) 



which, provided that the mean number of photons in the 
probe mode is small compared to iV, is also infinitesi- 
mal. This constraint is consistent with the assumption 
of negligible condensate depiction. 

Inserting these definitions into the equation of motion 
(E9) for CqCq, and keeping only terms linear in the in- 
finitesimal operators, gives 



l*"^"- 



(38) 



which has the trivial solution jg = 0. As a result, this op- 
erator can be dropped from the linearized equations for 
S- and S+ . This leads to a set of three coupled infinites- 
imal operators whose linearized equations of motion can 
be expressed as 



dr 



iM6, 



(39) 



where d = {6a, 6^,6+)'^ , the matrix M is given by 



M 



S ~x 


-X 


X (1 + /?) 


p 


'X ~(3 


-(1 + /3) 



(40) 



and we have introduced the dimensionless time r = u>rt, 
bJr = TiK^ /2ra being the atomic recoil frequency, as well 
as the dimensionless control parameters 



(41) 



2|A|w^ 



5 = 6' jlOr, 



and 



P 



AirfujNFo 



mujr 



(42) 



(43) 



|0) being the vacuum state. We proceed by linearizing 
the atomic field operators around their initial expecta- 
tion values, which can be determined from Eq. (pq), to- 
gether with the approximate commutation relations given 
by Eq. (|l|). This yields 



Here x is a dimensionless atom-probe coupling constant, 
5 is the pump-probe detuning in units of the atomic recoil 
frequency, and (5 gives the strength of collisions between 
the side modes. In Appendix B we give the effective 
Hamiltonian from which Eq. (B9) can be derived. 



The solution to Eq. (39) is then given by 



Sir) 



iMi 



^(0). 



(44) 



From this we see that the time dependence of the in- 
finitesimal operators is determined by the eigenvalue 
spectrum of the matrix M. For certain values of the 
control parameters Xi <^: ^iid /3, one of the eigenvalues 
contains a negative imaginary part. When this occurs, 
the infinitesimal operators undergo exponential growth. 
This exponential instability is the focus of the next sec- 
tion, where we discuss its properties in detail. 



VI. EXPONENTIAL INSTABILITY 



The eigenvalues of M are determined by the charac- 
teristic equation 

uj^ - Slo^ - (1 + 2(3)uj + (1 + 2(3)6 + 2x^ = 0, (45) 

which has either three real solutions, or one real and a 
pair of complex conjugate solutions. In the first case, 
the system is stable and exhibits only small oscillations 
around its initial state. In the second case, the system is 
unstable and grows exponentially, even from noise. 

From Eq. (^5|) one finds that exponential instability 
occurs when 



x'> 



[(3 + 60 + (52)3/2 + S^~- 95(1 + 2/3)]/27. (46) 



In Fig. 1 we plot the region of instability as a function of 
6 and x^- The shaded region of Fig.l corresponds to the 
instability region in the absence of collisions {(3 = 0). As 
collisions are added, the boundaries shift, illustrated by 
the dashed and dotted curves, which show the boundaries 
for the cases (3 = .3 and (3 = 1.0, respectively. 

From Fig. 1, we see that for positive 6 the boundary 
asymptotically reduces to 



X^ > 26^/27 -{1 + 2(3)6/6 



(47) 



i.e., the threshold value of x^ increases with the third 
power of the detuning and is only weakly influenced by 
the presence of collisions. On the other hand, for negative 
6 the asymptotic behavior is given by 



X^>-{l + 2p)6/2, 



(48) 



which only grows linearly with 6 and is strongly af- 
fected by interatomic collisions, which have the effect 
of reducing the unstable region. In earlier work [ p5[ , 
it was shown that this lower threshold occurs in the 
absence of collisions when atomic diffraction overcomes 
the bunching process. For positive scattering lengths. 




FIG. 1. Exponential instability region in the S-x^ plane. 
The shaded area gives the unstable domain in the absence 
of collisions (/3 = 0). The dashed and dotted curves show 
how the boundaries change as collisions are included. They 
correspond to the cases /3 = .3 (dashed) adn (3 = 1.0 (dotted). 



we note that formation of a density grating increases 
the mean field energy. Collisions, therefore, should join 
diffraction in opposing the bunching process, resulting in 
a higher threshold for the instability. For negative scat- 
tering lengths, on the other hand, bunching reduces the 
mean field energy, hence collisions should enhance the 
bunching process and oppose diffraction, thus lowering 
the threshold. Lastly, we note that the instability region 
touches x^ = when 6 = \/l + 2/3 « H- /?. This roughly 
corresponds to the conservation of energy in the scatter- 
ing of a pump photon into the probe by an atom intitially 
at rest. 

Once we have found the eigenvectors and eigenvalues 
of M, we can reexpress the solution (H) in the form 



6{t) = Ue^^^U^M(O) 



(49) 



where U is the matrix of eigenvectors of M, such that Uij 
is the zth component of the jth eigenvector, and Q is the 
diagonal matrix of eigenvalues, such that the ith diagonal 
element of il is the ith eigenvalue of M. In the unsta- 
ble regime, we have the eigenvalues cji, 102 = ^ + iF, 
and UJ3 = Q — iT, where wi and f2 are real, and F is 
real and positive. Thus lui corresponds to an oscillating 
solution, LU2 an exponentially decaying solution, and 073 
corresponds to an exponentially growing solution. Even- 
tually, this exponentially growing solution will dominate, 
at which time we can neglect the other two terms, yield- 
ing the approximate solution 



^.W-EOfe'5/c(0)e 



(r+iO)T 



(50) 



r(x,5,p=oo) 



r(x,5,p=0-3) 




r(x,5,p=i.o) 



FIG. 2. Exponential growth rate F as a function of the 
scaled pump-probe detuning S and coupling parameter x, for 
various values of the collision parameter /3. Figure 2a shows 
the case P — 0, while Figs. 2b and 2c show the cases P — 0.3 
and /3 = 1.0, respectively. 



where Qk — UjsU^f. . The range of validity for 
this approximation is roughly 1 < Ft <^ In(viV), 
where the lower limit is set by the requirement that 
the exponentially growing terms dominate, and the 




40 ^2 60 



FIG. 3. The dashed lines show the growth rate F as a func- 
tion of x^ for the case 5 — 1 + f3, with the values of /? specified 
in the figure. The solid line is the approximate expression 
given by Eq. (jsj) . 

upper hmit comes from the requirement that the side 
mode populations remain a small fraction of the total 
atom number. This condition therefore formaUy defines 
the exponential growth regime. 

The rate of exponential growth F has the explicit form 



2 
where 



and 



+ v? 



1/3 



-vV 



1/3 



-6(1 + 2/3) -x^ + 



27' 



q 



i's + ep + s^). 



(51) 



(52) 



(53) 



As this equation is complex and does not provide much 
insight, wc have plotted F as a function of 6 and x^ for 
three different values of /?. Figure 2a shows the limit of 
negligible collisions /? = 0, and Figs. 2b and 2c show the 
cases /3 = .3 and /? = 1 respectively. From these figures, 
we observe that the most significant effect of collisions 
is to shift the lower threshold. The values of F in the 
vicinity of the maximum (for fixed x), on the other hand, 
show less pronounced variations. 

In cases where x^ ^ I'^'^l, |/3| we have r « — x^ and 



V¥ 



+ r^ « X^- I^^ this case Eq. ( pl| ) reduces simply to 

F«V3(x/2)2/^ (54) 

Among other things, this shows that the gain scales as 
the number of atoms in the condensate to the 1/3 power. 
In Fig. 3, the growth rate F is plotted versus x^ with 
6 = I + (3, which roughly maximizes F for fixed x- The 
three dashed curves correspond to different values of the 
collision parameter /?, while the solid line gives the ap- 
proximate result (p4). This shows that the approxima- 
tion is a relatively accurate estimate of the maximum 
gain for all values of x^ ■ 



VII. QUANTUM STATISTICS 

In this section we use the solution ( p9| ) to compute 
some of the quantum statistical properties of the sys- 
tem. This, however, first requires a more detailed dis- 
cussion of the physical meaning of the infinitesimal op- 
erators. The first, Sa = a/vN, is clearly just a rescaling 
of the photon annihilation operator. From it one can 
compute all properties of the electric field and/or the 
photon satisfies of the probe mode. The atomic side 
mode operators S- — c_cq/N and ^+ = CqC+/N, how- 
ever, are not simply rescalings of atom annihilation op- 
erators. Rather, they are directly related to the atomic 
density, p{r) = ^l{r)^g{r). 

To illustrate this point, we expand p(r) according to 
Eq. (ETI) and linearize, yielding 



p(r)-iV|^o(r)| = 



.K.r (^ 



6^ 



H.c. 



(55) 



From this expression we see that the side mode operators 
S± indeed describe the appearance of a density modula- 
tion with wavelength 2'!r/K . 

In addition to the atomic density, one can also express 
the number operators for the side modes in terms of S± . 
For example, we have after linearization 



NS-Sl = 



N 



4, jN+l) 
c c_ . 



(56) 



Hence with N + 1 ^ N the number operator for the '-' 
side mode can be expressed as 



cLc- 



nLsI 



Similarly the number operator for the '- 
given by 



(57) 



side mode is 



,lc 



N6U 



+ ■ 



(58) 



From these number operators, one can therefore com- 
pute the number statistics of the side modes in the linear 
regime. 

From the analytical solution (E^) it is straightforward 
to compute the properties of the atomic and optical fields 
for an arbitrary initial condition. We focus on two condi- 
tions which appear readily accessible experimentally. In 
the first one, the probe field and the atomic side modes 
all begin in the vacuum state. In this case the exponen- 
tial growth is triggered by vacuum fluctuations in both 
the probe field and the atomic density. A second pos- 
sible triggering mechanism involves injecting of a weak 
laser field into the probe mode. Both initial situations 
are investigated by assuming that the probe mode is ini- 
tially in the coherent state a, such that a|a) = a\a), the 
vacuum case corresponding to a = 0. In addition, we 
assume throughout that the condensate side modes be- 
gin in the vacuum state. Hence, the initial state of the 



three-mode system can be expressed as |a, 0,0), where 
the first index refers to the probe mode, and the second 
and third indices give the states of the momentum side 
modes. 



A. Electric field and atomic density 

The expectation value of the operator 5a is sufficient to 
compute the mean electric field, and likewise, the mean 
values of (S_ and 5+ are sufficient to compute the mean 
atomic density. We now give analytic solutions for these 
physical quantities and their quantum mechanical uncer- 
tainties in the exponential growth regime, where all but 
the leading exponential terms can be safely neglected. 

The electric field operator for the probe field is given 

by 



E{r) = e{r)£ir,T)VNSair)+H.c., 
where e'(r) is the polarization unit vector, and 



f(r,T)=-i 



2^|f^o||A| 



(pB(r)e~'("''/'^'-)^ 



(59) 



(60) 



contains all constants of proportionality, the normalized 
spatial wavefunction of the probe mode (/'^(r), and the 
oscillation at the pump frequency ujq. The mean electric 
field is obtained by inserting Eq. (pQ), and taking the 
quantum mechanical expectation value with respect to 
the initial state |q;,0,0), which yields 



e{r) ■ (E{r)) = £{r, rXaaCce^^^'''^^ + c. 



(61) 



This corresponds to an oscillating mean field with ampli- 
tude 



ii;o(r) = 2|£(r,T)||Caa||a|e 



Ft 



(62) 



From this expression, we see that there is a nonzero mean 
field provided only that a 7^ O. We also see from Eq. ( |6f| ) 
that the mean field amplitude grows exponentially at the 
rate F, and that its frequency is shifted by — f2 away 
from the pump frequency. Its phase, on the other hand, 
has a somewhat complicated dependence on the system 
parameters. An analytic expression for this phase can 
be computed directly from Eq. (|6l[), but as we draw no 
specific conclusions from it, we do not give the explicit 
expression here. 

The variance in the electric field can also be computed 
in a straightforward manner from Eq. (pO|), yielding 



Aii;(r) = \/2|f(r,T)||Ca-|e 



Tt 



(63) 



This shows that the fluctuations also grow exponentially 
in time, irrespective of whether the mean field vanishes 
or not, and are in fact independent of a. Hence these 
fluctuations can be attributed solely to the amplification 
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of quantum noise, i.e. vacuum fluctuations in the probe 
electric field as well as atomic density fluctuations. While 
the mean field and the fluctuations both grow exponen- 
tially in time, the relative uncertainty, on the other hand, 
is constant in time, given by 



A£;(r) f{S,x,f3) 



Eoir) 



V2\c 



(64) 



Here, we have introduced the fluctuation function 



f{S,X,f3) = 



\u,--\ 



(65) 



which can be computed directly from the eigenvectors of 
the matrix M, and therefore depends only on the pa- 
rameters X: <5, and (3. Figure 4 plots f{S, x, (3) above the 
S-x^ plane for various values of the collision parameter (3. 
Figure 4a shows the limit of negligible collisions {/3 — 0), 
where we see that f{S,Xj 0) is nearly flat in the vicinity of 
maximum gain [S ~ 1), where it has a value somewhere 
between 1 and 2. It steadily increases from this value 
as the pump-probe detuning S moves in the negative di- 
rection. Figures 4b and 4c show the cases (3 — 0.3 and 
/3 = 1.0 respectively. From these we see that the effect of 
increasing the collision parameter is to flatten f{S,x,f3) 
as a function of x and 6. 

In a similar manner, we next calculate the mean value 
and variance of the atomic density. By inserting Eq. 
( pO| ) into Eq. (psh and taking the expectation value with 
respect to the initial state, we find the expectation value 
of the atomic density to be 

(p(r)) = iV|(po(r)P + Po(r, r) cos[K • r + fir + 0] (66) 



where the amplitude Pq(v) of the density modulation is 
given by 



Po{r,T)^VN\Mr)r\C-a + C+a\\a\e 



Ft 



(67) 



Thus we see that the mean atomic density is the sum of 
two contributions, the initial density of the condensate 
plus a density modulation which grows exponentially in 
time, provided of course that a ^ 0. Together with Eq. 
( pi] ) , this shows that the phase symmetry of the system 
is broken by the phase of the injected field. Only in 
the case a = does the symmetry remain unbroken. We 
note that for the atomic side modes, this is not symmetry 
breaking in the commonly used sense of non-zero mean 
fields. Rather, it is the mean atomic density modulation 
which acquires a nonzero phase. Note also that the mean 
density modulation is not stationary, as its phase is given 



f(X,5,|3=00) 



f(X,5,|3=0-3) 




f(X,5,|3=1.0) 



FIG. 4. Comparison of the the fluctuation function 
f{5, Xi P) as a function of 5 and y^ for various values of the 
collision parameter 13. Figure 3a gives the limit of negligi- 
ble collisions (/3 = 0), where Figs. 3b and 3c show the cases 
/3 = 0.3 and j3 — 1.0, respectively. 

by Qt + (j), where both il, and (f) are depend only on the 
system parameters 6, x, and f3. 

The variance in the atomic density can also be read- 
ily computed in the exponential growth approximation, 
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yielding 



Apir) = 2VN\Mrr\C— + C+-\e' 



(68) 



which shows that the density fluctuations grow exponen- 
tiaUy in time, even in the case a = 0. The ratio between 
the density variance and the modulation ampHtude is 
constant in time, and is given by the same expression as 
that for the probe, i.e.. 



Ap(r) f{S,x,P) 



Po(r) 



V2\c 



(69) 



From Eqs. ( |64| ) and (|69|), we see that in the case 
f{5,x,l3) <C |a|, both the mean electric field and the 
mean atomic density modulation arc quantum mechan- 
ically well-defined, meaning that the quantum noise is 
small compared to their mean values. In this regime, 
both quantities could be adequately treated as classical 
(c-number) fields. Outside of this regime, however, the 
quantum fluctuations play a significant role, and a clas- 
sical description no longer suffices. 

The main implication of these results is that by vary- 
ing the system control parameters, and in particular the 
injected field intensity and phase, one can vary the mean 
electric field and atomic density modulation continuously 
between two limits. For |a| <C f{S,X:P) the fields are 
dominated by quantum fluctuations, and we can expect 
to find important non-classical effects. In the limit of a 
'strong' injected field, however, the fluctuations are not 
signiflcant, and the atomic and optical flelds behave clas- 
sically. 



B. Intensities 



We now turn to the number statistics of the three field 
modes, concentrating on the mean atom/photon numbers 
and their variances. It is convenient to r eexp ress the 
mode number operators given by Eqs. (|56|)-(|5q), as 



N, 



NS]6, 



5. 



(70) 



where the index j is again the mode label a, — , or -f , 
and the (5-function accounts for the fact that we have 
normally ordered the infinitesimal field operators. It is 
then straightforward to derive the full time-dependent 
solution for the mean occupation numbers Nj = (Nj) as 



N, = \a,aiT)\M^ + \aj.{T)\^-Sj 



where 



,M-E^' 



^kU^e^'^'^- 



(71) 



(72) 



ln[N, 




FIG. 5. Logarithmic plot of the probe intensity Na as a 
function of time. This thick curves show the exact solution, 
given by Eq. (iTlJ), and the corresponding thin lines give the 
approximate solution of Eq. (uS). The parameters chosen are 
S — 1, x^ ~ 1] S'lid (3 — 0. Each pair of curves corresponds 
to a different value of of the initial probe intensity \af , as 
specified in the figure. 

The first term in Eq. ( [7l|) can be interpreted as the 
stimulated contribution to the intensity. 



WjU = Wja{T)\^\a\' 



(73) 



whereas the second term gives the spontaneous contribu- 
tion. 



m,,-\a,-(r)\' 



(74) 



present even in the case a = 0. In the exponential growth 
regime Eq. ( [7l| ) reduces to 

(75) 



N,^{\Qa\'\a\' + \Q-ne'^\ 



which shows that the mode occupation grows exponen- 
tially, even in the spontaneous case, where it was seen 
that the mean electric field and mean density modula- 
tion both vanish. The validity of the exponential ap- 
proximation ( [75|) is demonstrated in Fig. 5, where we 
have plotted the logarithm of the probe intensity as a 
function of time. The parameters chosen for the plot are 
S = 1, x^ ~ 1, and (3 = 0. The thick lines give the full 
solution (^ ) for three different values of the initial probe 



intensity \a 



0, 1, 10. The corresponding thin lines are 



the approximate solutions given by Eq. (^q), which sow 
good agreement for r > 1. 

Turning now to quantum fluctuations in the occupa- 
tion numbers, we find that the relative uncertainties are 
given by 






= Jl- 



\aja{T)a\ 



1 



(76) 



fe=i 



From Eq. (WWj it is clear that the second term un- 
der the radicand in the above expression is < 1, 
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FIG. 6. The side mode number variance AN^/N^ is plot- 
ted versus time (thick solid line) . Also shown are the approx- 
imate solution given by Eq. ( |77| ) (thick dashed line) as well 
as the variances for a thermal state (thin solid line) and a 
coherent state (thin dashed line) with the same mean value 
A''-. The parameters chosen are 5 = 1, x^ = i, (^ ~ 0, and 




FIG. 7. The long time limit of ANj/Nj is plotted against 
the initial probe intensity |a|^ (thick solid line). It various 
continuously between thermal (upper dashed line) and coher- 
ent (lower dashed line) limits. The parameters chosen are 
A = 1, x^ = 1, and /3 = 0. 

which means that as a function of time, the relative un- 
certainty is always between y^l/N, characteristic of the 
fluctuations found in a coherent state, and y^l + 1/N, 
which is the signature of thermal number fluctuations. 

While for very short times the relative uncertainty may 
fluctuate between the thermal and coherent limits, once 
the exponentially growing terms dominate, the relative 
uncertainty eventually reaches a steady-state value given 
by 






Thus we see that when |ap :$> /^(x, 5, /3), the relative un- 
certainty tends towards zero, while in the opposite case, 
it tends toward one. These limits can be labeled as the 
stimulated and spontaneous limits respectively. In the 
intermediate regime, the fluctuations can be varied con- 
tinuously between the thermal and coherent limits, e.g. 
by varying the injected laser intensity, thus achieving op- 
tical control over the quantum statistics of matter waves. 
The behavior of the particle number variances is il- 
lustrated by Figs. 6 and 7. In Fig. 6 the full time 
dependence of AN-/N- is shown (thick solid line). Also 
shown are the approximate solution given by Eq. ( [77| ) 
(thick dashed line) as well as the variances for a thermal 
state (thin solid line) and a coherent state (thin dashed 
line) with the same mean value N- . The parameters cho- 
sen are J = 1, x^ = Ij /? = 0; and a = 1. Thus we see 
that the variance always falls between those of thermal 
and coherent fields. We also see that the long time be- 
havior is well approximated by Eq. (^). In Fig. 7 the 
steady state value for large r is plotted as a function of 
the initial probe intensity |ap. Thus we see that it is 
possible to vary the output continuously over the whole 
range between thermal and coherent limits, simply by 
varying the initial probe intensity. The parameters cho- 
sen for the figure are 5 — 1, x^ — ^, and (3 — 0. 



VIII. ATOM-PHOTON ENTANGLEMENT 

We have previously discussed the analogy between the 
present system and the non-degenerate optical paramet- 
ric amplifier (OPA). One of the most interesting applica- 
tions of the OPA is the generation of entangled quantum 
optical states. We show that similar entanglements occur 
in the present system, but they are now between atomic 
and optical field modes. We first examine the two-mode 
intensity correlation functions, which give a measure of 
entanglement, and can be used to determine whether or 
not non-classical correlations exist between the three field 
modes. We then discuss the issue of two-mode squeez- 
ing, and show how this phenomenon manifests itself in 
the present system. 



A. Two-mode intensity correlations 

The equal-time intensity correlation functions are de- 
fined in the usual manner as 



(2) 

9ij 



mm 



(78) 



For classical fields, the two-mode (i ^ j) correlations are 
constrained by the Cauchy-Schwartz inequality 
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FIG. 8. The cross-correlation function g_^ is plotted as 
a function of time (solid line). The upper dashed line gives 
the maximum allowed by quantum mechanics, while the lower 
dashed line gives the upper limit for classical correlations only. 
The parameters chosen are 5 = 1, x^ = 1, and /3 = 0. Figure 
8a corresponds to triggering from noise (a — 0), while Fig. 
8b shows the case a = 1. 

Quantum mechanical fields, however, can violate this in- 
equality and are instead constrained by 



(2) ^ 
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which reduces to the classical result in the limit of large 
intensities. 

We focus our attention on the spontaneous case a = 0. 
Here the single-mode intensity correlation functions are 
those of thermal fields, gl (r) = 2. In this case, the 
equal-time intensity cross-correlation functions are found 
to be 
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violate the Cauchy-Schwartz inequality, while g^Hr) 
is consistent with classical cross-correlations. Further- 
more, the explicit evaluation of the Cii's shows that 

I+{t) <C Ia{T), which implies that gali^) is very close 
to the maximum violation of the classical inequality con- 
sistent with quantum mechanics, whereas for g_^(T) the 
violation is not close to the allowed maximum. In the 
two-mode parametric amplifier, the two-mode correlation 
function shows the maximum violation of the Cauchy- 
Schwartz inequality consistent with quantum mechanics. 
In the three-mode system, however, the two-mode cross- 
correlation functions involve a trace over the third mode, 
hence it is not surprising that they are not maximized. 
This is illustrated in Fig. 8a, where we have plotted the 
correlation function g_^ as a function of time (solid line). 
Also shown for comparison are the quantum mechanical 
upper limit given by Eq. (pG) (upper dashed line), as 
well as the classical upper limit given by Eq. ( [79[ ) (lower 
dashed Hue). The parameters are 5 — \, x^ — 1, (3 — Q, 
and a — Q. We see that when the system is triggered by 

(2) 

vacuum fluctuations alone, g_a is virtually indistinguish- 
able from the maximum allowed by quantum mechanics. 
If we now allow for an injected coherent probe field 
{a 7^ 0), we first note that the intensities are increased 
by approximately jap, which means that the time scale 
on which the classical and quantum upper limits ( [79[ ) and 
( pO| ) converge is reduced by l/|ap, making an experimen- 
tal confirmation of quantum correlations more difficult. 
In addition, whereas for the spontaneous case a — Q, 

(2) 

numerics show the cross-correlation g^_ follows almost 
exactly the quantum upper limit (pG) for all t > 0, for 
a 7^ 0, it lies somewhere in between the quantum ( |80| ) 
and classical (79) limits. As a is increased, it falls ever 
closer to the classical upper limit, so that in the limit of 
very large a, the fields exhibit classical cross-correlations 
only. This effect is illustrated by Fig. 8b, which is iden- 
tical to 8a, except that now we have taken a = 1. This 
shows that even when the probe field initially contains 
only one photon on average, it is enough to significantly 
reduce the quantum correlations from the quantum me- 
chanical maximum. 



B. Two- mode squeezing 

To complete the study of atom-photon entanglement 
we complement the investigation of intensity cross- 
correlations with a discussion of phase-sensitive two- 
mode correlations. Drawing again on the similarity be- 
tween the current model and the two-mode parametric 
amplifier we expect the correlations between the cavity 
mode and the atomic side-mode "— " to be of particu- 
lar interest as a "squeezing-like" behavior may occur. In 
analogy to the parametric amplifier we thus introduce the 
quadrature components of the superposed atomic and op- 
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tical fields 
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h.c. 



(82) 



The variance of Xg is given by 

V{Xg) = {XD - {Xgf 

^\aa-{T)e'''+a^^{T)e-''\\ 



(83) 



Note that this variance ( |83| ) is independent of the injected 
signal strength a, just as with the electric field, atomic 
density, and mode intensities. It follows from Eq. ( [8^ ) 
that the angle 9min{T) that minimizes the quadrature 
variance is determined by 



arg[aa-(T)al_(r)] + 29„iiniT) = tt 



(84) 



so that the corresponding minimum of V for fixed r is 
given by 



V;n»(T) = [|aa-(T)|-|a— (r)|] 



(85) 



The Heisenberg uncertainty principle gives 
V {Xg)V {Xgj^Tj /2) > 1, hence, a quadrature component 
is squeezed provided V{Xo) < 1. 

In Fig. 9, we display VminiT) for various values of the 
system parameters. We find that for 5 « 1, i.e., maxi- 
mum exponential growth rate, Vmin is a concave function 
of T displaying a single (global) minimum which is typi- 
cally of the order of 10^^ (cf. the full and dashed curves). 
The maximum squeezing time r^, i.e., the largest t for 
which VminiT) = 1, is given by t^ = 3.5 for x = 1, ^ = 1, 
and (3 — whereas for x = 10 Tm decreases to a value 
of 0.20. The reduction is due to the increase in the ex- 
ponential growth rate T with larger x- From Fig. 9, we 
see that squeezing does indeed occur over a broad range 
of parameters, however it only persists over intermedi- 
ate time scales. For long times, Vmin{T) is dominated by 
the exponential behavior, which eventually leads to the 
violation of the squeezing condition. This is in contrast 
to the two-mode OPA, where the quadrature component 
remains squeezed for all time. 

To understand how the presence of a third mode 
quenches the squeezing, we reexpress the squeezing con- 
dition with the help of Eq. ([74]), yielding 



([7V„],p-[iV^]j'<2[7Vj 



(86) 



As our analysis of squeezing so far has not made ex- 
plicit use of the particular form of the matrix M, 
it can also be applied to the standard two-mode 
optical parametric amplifier. One simply chooses 
a suitable M where the third mode is decoupled. 
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FIG. 9. Minimum variance Vmin. of the quadrature com- 
ponents Xe as a function of r for parameter values x = 1-0, 
5 = 1.0, /3 = 0.0 (full curve); x = 3.0, S = 1.0, /3 = 0.0 
(dashed); x = 3.0, S = -17.0, P = 0.0 (dotted). 

In this case, the solution is well known [Q, and the two 
modes are symmetrical, with equal populations. This 
means that the l.h.s. of Eq. ( ^q) is always zero, and the 
squeezing condition is satisfied for all times. However, 
once the third mode is included, it introduces a small 
imbalance between the populations of the two original 
modes, since momentum conservation in requires that 
Na + N+ = N-. On long time scales, the l.h.s. of Eq. 
( pq ) grows like exp(i4rT), whereas the r.h.s. only grows 
as exp(i2rr). Therefore, the squeezing condition (Isl) 
must eventually cease to be satisfied. 

The introduction of moderate coUisional interactions, 
i.e., g < 1, only leads to quantitative modifications of 
Vmin without altering the characteristic features. An in- 
teresting change in behavior occurs, however, if S is tuned 
closer to the borders of the amplification range (dotted 
curve in Fig. 6). Under these circumstances squeezing 
still occurs but Vmin displays an oscillating behavior as 
a function of r. These oscillations are caused in part by 
the reduction of the exponential growth rate and a simul- 
taneous increase of the imaginary part of the eigenvalues 
of the matrix M. 

The extremal angle Omin varies as time evolves, even- 
tually attaining a constant value when the behavior of 
Vmin{T) is dominated by exponential growth. However, 
in many cases V{Xg) is well approximated by Vmin{T) if 
9 is chosen in the vicinity of 9min (t) ■ 



IX. DISCUSSION AND CONCLUSION 

The fundamental time scale in the system is the in- 
verse growth rate F"^. It is estimated by Eq. ( [5^ ) 
in units of the inverse recoil frequency (jJ~^ which, for 
sodium is given by LOr «1.7 fis. The estimated growth 
rate is then of the order of x^^^ — iV^n^g/ujr)'^^'^, 
where Ue = (/o/8/sat)(7/^)^- Here Iq is the pump in- 
tensity, Isat the atomic saturation intensity, and 7 the 
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atomic spontaneous decay rate. For sodium atoms we 
have Isat = 6.33mW/cm^. The parameter Uf, equals the 
fraction of excited atoms, hence, under the far-off reso- 
nance conditions we are considering, we have rig ^ 1. 
If we chose the ring cavity length L ~ 0.1m and effec- 
tive cross section S = lO^^m^, the atom-cavity coupling 
constant g is of the order 10^ s^^ so that g/ujr « 1. As 
one has a great latitude in choosing the values of rie and 
the total number N of atoms in the BEC it should be 
possible to vary the exponential growth rate over a wide 
range. 

Our choice of n and Ug is constrained, however, by 
the requirement that the spontaneous heating rate C be 
much smaller than the exponential growth rate F, so that 
spontaneous emission can in fact be neglected. The heat- 
ing rate (in units of tOr) is given by £ = Uej/uJr- As 
j/ujr ~ 10 the condition F ^ £ translates into N ^ 10'^ 
n'^. As Tie <C 1, this condition is practically always ful- 
filled. Also related to spontaneous emission is the two- 
body dipole-dipole interaction, which acts in addition to 
ground-state collisions. For very cold atoms whose de 
Broglie wave length is large in comparison to the pump 
laser wavelength the dipole-dipole interaction can be ap- 
proximated as a contact interaction pS] . Comparing the 
strength of this interaction to that ground-state collisions 
one finds that the former is negligible under the condi- 
tion He ^ SfikQa/m^, which translates to rie ^ 10~^ for 
sodium atoms. If this condition is not met the dipole- 
dipole interaction can still be accounted for to a good de- 
gree of approximation by modifying the scattering length 
according to a — > a — IqJ^ /32Isat^'^(^r- 

Another important parameter is the collision parame- 
ter /3. Estimating the quantity NFq to be of the order 
of the atomic density we find the collision parameter /? 
to lie in the range 0.1—1 for 'typical' densities around 
10^^ cm^^. If /3 could be measured by, e.g. observing 
the boundaries of the instability regime, than this could 
be used as a novel means to determine the atomic s-wave 
scattering length. 

In a realistic optical cavity, the lifetime of a photon is 
of the order of 1-lOns, corresponding to a decay rate of 
10^ to 10"* in units of the recoil frequency. This tells us 
how large the growth rate F would have to be to result in 
a buildup of photons in the cavity. Since F goes like N^^'^ 
we would likely need a very large condensate, with 10^^ or 
more atoms, to achieve this. In the future, larger conden- 
sates and better cavities, will be available, at which point 
the theory could presumably be tested. Future research 
will study the quantum statistics of the system under the 
influence of probe field damping, which will relate the 
model more closely to current experiments. For exam- 
ple, recent experiments by W. Ketterle's group involving 
a condensate driven by a far-off resonant pump laser have 
demonstrated the apperance of distinct momentum side 
modes as a consequence of spontaneous emission. While 
these experiments involve at most one spontaneous pho- 
ton at a time in the condensate, they clearly demon- 
strate many aspects of the theory we have described, as 



wll as the importance of continuing to develop a nonlin- 
ear/quantum optics approach to the theory of optically 
driven Bose-Einstein condensates. 
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APPENDIX A: NON-ORTHOGONALITY 

By properly taking into account the non-orthogonality 
of the atomic field modes, it can be shown that the only 
surviving effect in the linearized theory is the modifi- 
cation of the atomic polarization term in the equation 
of motion for the probe field ( [lO|) to include a second 
scattering mechanism in which a condensate scatters a 
photon without changing its center of mass state. As 
a consequence of momentum conservation, this process 
is supressed by a factor {(po\ip-) relative to the process 
which transfers the atom to the side mode state. Bose 
enhancement, on the other hand, is stronger for this tran- 
sition by a factor vN, because we now have N identical 
bosons in both the initial and final states. Thus it is the 
product VN{(po\ip-) which must be negligible if we are 
to make the orthogonality approximation. 

More precisely, we find that the quantum statistics 
must be modified by taking a — > a — s, where 



(fi - iV) 



(87) 



in order to account for this additional scattering mech- 
anism. Thus when a is comparable to /(i5, x,/3), the 
condition \s\ <C |a| allows us to neglect s. For the 
case a <C /((5, x,/3), one the other hand, requires \s\ <ti 
/((5, x,/3). Since f{S,x-f3) is typically of order 1, then the 
condition \s\ <^ 1 is sufficient to satisfy both conditions. 
For a dilute condensate, the ground state wave function 
<y9o(r) is given to a approximately by the single-particle 
ground state. For a harmonic trap with K taken along 
the z-axis, one has thus 

(^o|^->=e-I(^^^)^ (88) 

we K approximately twice the optical wave number, and 
Wc is the condensate width along K. For a dilute con- 
densate which is an order of magnitude or larger than 
the optical wavelength, this integral is clearly vanishingly 
small, and s « for all reasonable values A^ and the sys- 
tem parameters (5, x, and /3. 

In the case of a dense condensate the density distribu- 
tion can be described with the help of the Thomas-Fermi 
approximation, i.e.. 
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|0o(r)P - [m - K,(r)]m/(47r?iViV), 



(89) 



where the chemical potential /i is determined from the 
normalization requirement. Starting from Eq. (|89| ) and 
choosing the recoil wave vector K to lie parallel the pr- 
axis, the overlap integral is found to be 



((^olv'-) = 1532{KW,)/{KW,), 



(90) 



where J2 is the modified Bessel functio, and the W is 
the condensate width along K. For condensates which 
are large compared to the optical wavelength, Eq. ( |90| ) 
behaves like —15 sm{KWc) / (KWc)^ . So assuming that 
KWc = 100, which corresponds to a lO/im condensate, 
we would need x-\/iV/(ri — iV) to be of the order 10^ 
in order achieve an appreciable value of s, and this is 
much larger than is currently feasible. However, if, for 
example, K is oriented parallel to a trap axis with tight 
confinement KWc need not be too large and diffraction 
effects may be appreciable. 

A more eflticient way to increase s consists in confin- 
ing the condensate in three-dimensional rectangular trap 
potential. In this case the Thomas-Fermi approximation 
gives |(/9o(r)P = 1/V^, V being the volume of the trap. 
With K chosen along the z-axis of the trap, and Wc be- 
ing the width of the box in this direction, one obtains for 
the overlap integral 



tKWj2 ^HKWc/2) 

KWc/2 



For large KWc I now decays only as {KWc)~^. 
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APPENDIX B: EFFECTIVE HAMILTONIAN 

If we define C- = \/N6_ , and (7+ = \/N6+ , then these 
operators obey approximately bosonic commutation re- 
lations, and the system of equations (|3^) can be derived 
from the effective Hamiltonian 

Heff = (1 + l3)ClC- + (1 + P)ClC+ - a^a 



+ (3{Clcl + C-C+) 

+ xia^'Cl + aC- + a^C+ + aC|). 
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